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Recently BCS superconductivity at 203 K has been discovery in a highly compressed hydrogen sulfide. We 
use first-principles calculations to systematically examine the effects of partially substituting the chalcogenide 
atoms on the superconductivity of hydrogen chalcogenides under high pressures. We find detailed trends of how 
the critical temperature changes with increasing the V-, VI- or Vll-substitution rate, which highlight the key 
roles played by low atomic mass and by strong covalent metallicity. In particular, a possible record high critical 
temperature of 280 K is predicted in a stable H3S0.925P0.075 with the ImSm structure under 250 GPa. 


I. INTRODUCTION 


Since the discovery of superconductivity in mercury at 4 K 
in 1911', scientists have been ardently pursuing new and bet¬ 
ter superconductors at higher temperatures. In a conventional 
superconductor, the vibrations of crystal lattice provide an 
attractive force that binds an electron with its time-reversal 
partner into a Cooper pair’. The Cooper pairs can Bose con¬ 
dense below a critical temperature (Tc), which had been be¬ 
lieved to be no larger than 25 K prior to the discovery of 
MgB 2 ^’^. In this case, the Debye temperature could be as large 
as ~ 10^ K, and low Tc is primarily related to small electron- 
phonon couplings^. Since 1986 and 2008, respectively, the 
discoveries of copper- and iron-based superconductors have 
provided two new avenues for making high-T^ superconduc¬ 
tors^*^, while generating new excitements in fundamental 
physics. Although their unconventional mechanisms are still 
under hot debate, to date a record Tc of 164 K has been exper¬ 
imentally realized in the cuprate family® and 56 K among the 
iron-pnictide compounds*^. 

One may thus wonder whether a Tc higher than 164 K 
could be achieved and whether the conventional phonon- 
mediated mechanism could play a significant role in such a 
race. Both answers have become extremely positive and fasci¬ 
nating, given the recent theoretical prediction*^’*'* and partic¬ 
ularly the most recent experimental observation*^’*® of super¬ 
conductivity at 203 K in a highly compressed hydrogen sul- 
fide*^“^^. As pointed out by Bernstein et al. in a microscopic 
theory of H 3 S*®, the substantial covalent metallicity leads to 
a large electron-phonon coupling, and the low atomic mass 
leads to high-frequency phonon modes. Both features play es¬ 
sential roles in increasing the Tc of H 3 S. Notably, the former 
feature is analogous to that of MgB 2 ^®’^® whereas the latter 
effect is similar to those in the H-rich materials^®^^^’^®’ 

As a powerful method for the optimization of Tc, the 
atomic substitution has been routine in superconductivity ex¬ 
periments^®’'***’'*^’ . However, there is yet any study on how 
the atomic substitution influences the Tc of H 3 S and of the 
more general hydrogen chalcogenides. Here, we systemati¬ 
cally examine the effects of (partially) substituting the chalco¬ 
genide atoms on the superconductivity and particularly the 
Tc’s of hydrogen chalcogenides in the /m3m phase under 
high pressures, based on the first-principles calculations with 
virtual crystal approximation (VCA)'*-*^®. In the V- and VII- 


substitution cases, we find that the significant changes of the 
DOS at Fermi surface and of the phonon linewidths, coming 
from the different number of valance electrons, are the princi¬ 
pal factors affecting the electron-phonon coupling and Tc- In 
the particular case of H 3 Si_xPx at 200 GPa, we find that the 
DOS, the electron-phonon coupling constant, and the Tc all 
first increase and then decrease as the P-substitution rate in¬ 
creases from zero to 0.2. In the optimized condition in which 
X = 0.075 and the pressure increases to 250 GPa, a possi¬ 
ble record high Tc of 280 K is predicted. In contrast, the Tc 
decreases monotonously with increasing the Vll-substitution 
degree. In the Vl-substitution cases, the Tc does not appear to 
increase substantially, because of the reduction of DOS with 
O- or Se-substitution, or due to the softening of logarithmi¬ 
cally averaged phonon frequency with Te-substitution. These 
findings emphasize the importance of low atomic mass and 
strong covalent metallicity in conventional high-Tc supercon¬ 
ductors, pave the way for substantially enhancing Tc by com¬ 
bining application of a high pressure and properly designed 
chemical substitution, and suggest that in principle it is not 
impossible to boost Tc to ice point in optimized conditions. 


II. RESULTS 

Soon after the experimental report*®’*® on the record Tc of 
a highly-compressed hydrogen sulfide*®’*'*, further theoreti¬ 
cal studies*®^®® have reached two consensuses: (i) the super¬ 
conducting matter is most likely H 3 S with a strong covalent 
character in the Im3m phase, as sketched in Fig. la, and (ii) 
the superconductivity of H 3 S is phonon-mediated. Thus, the 
superconductivity of highly-compressed hydrogen sulfide can 
be accurately described by the celebrated Eliashberg theory'*®, 
which takes into account the renormalization of electron- 
electron repulsion by electron-phonon interactions. In this 
theory, the Allen-Dynes-modified McMillan formula'*®’®** re¬ 
lates Tc to the logarithmically averaged phonon frequency win, 
the effective Coulomb repulsion pf, and the electron-phonon 
coupling constant A: 


Tc = /1/2 


Wln_ 

1.20 


exp 


1.04(1 +A) 

X^ 7 ?(T*+a 62 Ay 


( 1 ) 


where /i and /2 are the strong coupling and the shape cor¬ 
rection factors'*®, respectively. Clearly, there are two ways to 
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FIG. 1. ImSm structure of H 3 S and phonon spectra with phonon linewidth of Vl-suhstitution systems, (a) The ImSm structure of 
hydrogen chalcogenides. (b)-(e) Phonon spectra and phonon linewidths of four representative Vl-substitution systems: (b) H3S, (c) HsSe, 
(d) H3S0.4O0.6, and (e) £[380.7X60.3. The magnitude of the phonon linewidth is indicated by the size of red error bar, and the magnitude for 
£[ 380 , 400,6 is plotted with twice of the real values. 


enhance T^, i.e., to search for materials with high-frequency 
phonon modes and to increase electron-phonon couplings. As 
mentioned in the introduction, the H 3 S happens to consti¬ 
tute both advantages and hence has a record high T^. In¬ 
deed, within this framework recent theoretical calculations 
have well explain the high Tc in the highly-compressed hy¬ 
drogen sulfide*^”^^. Therefore, we will first reproduce the Tc 
of H 3 S in the ImZm phase at 200 GPa and study the substitu¬ 
tion of S atoms by other VI atoms. Comparing our results with 
previous ones will show the validity of our first-principles cal¬ 
culations with VC A. 

In the Vl-substitution systems at 200 GPa, we focus on 
studying H 3 Si_xSex in the full range of x, H 3 Si_xOx in the 
range of a; = 0.0 ~ 0.6, and H 3 Si_xTex in the range of 
X = 0.0 ^ 0.3. Because when x takes a greater value in the 
case of O- or Te-substitution, we can find imaginary phonon 
modes at T point, indicating a lattice structure instability. The 
electronic band structures of these systems are shown in the 
supplementary materials, and our results for H 3 S and HsSe 
are in great agreement with previous reports based on sim¬ 
ilar first-principles calculationsWe find that the VI- 
substitutions have little influence on the electronic band struc¬ 
ture around the Fermi energy, except unimportant changes 
near F point (see supplementary materials). We further And 
that including spin-orbit couplings for the Te-substitution case 
does not introduce any significant correction, either. Notably, 
the main electronic effect is the reduction of DOS at Fermi 
surface in all of these Vl-substitution cases, as summarized 
in Table I. The DOS decreases monotonically with increasing 
the Vl-substitution rate. Moreover, the lighter the substitution 
element, the stronger the reduction of DOS. 

Figures Ib-le display the phonon spectra of four repre¬ 
sentative Vl-substitution systems, and the magnitude of the 
phonon linewidth is indicated by the size of the red error bar. 
Our results for H 3 S and £[386 are again consistent with the 


earlier reports Overall there is a clear separation between 
H modes at high energy and S modes at low energy*®. Ev¬ 
idently, with increasing the degree of heavy-element (Se or 
Te) substitution, all the three acoustic phonon modes decrease 
in frequency. In addition, the H-VI bond-bending modes (dis¬ 
placement of an H atom perpendicular to a H-VI bond) in the 
intermediate frequency region are softened, whereas the H-VI 
bond-stretching modes (displacement of an H atom parallel 
to a H-VI bond) in the high frequency region are hardened. 
Not surprisingly, the opposite trends occur for the case of O- 
substitution. The above changes at low and intermediate fre¬ 
quency are expected by considering the relative atomic mass 
of the substitution elements, whereas the changes at high fre¬ 
quency might be explained by the dependence of the chem- 


TABLE I. Fermi-surface DOS (in units of Hartree“*^/spin) of hy¬ 
drogen chalcogenides in various V-, VI-, and Vll-suhstitution 
systems at 200 GPa. 


X 

H3Sl a:02: 

H 3 S 1 —xSCx 

H3Sl-a:Te2, 


0.0 

7.91 

7.91 

7.91 


0.1 

7.70 

7.78 

7.90 


0.3 

6.84 

7.61 

7.53 


0.6 

4.53 

7.51 



1.0 


7.31 



X 

H3Sl a;Pa; 

H3S6i x 

H3Si,CE 

H3Sei 3;;Brx 

0.000 

7.91 

7.31 

7.91 

7.31 

0.025 

8.57 

7.80 

7.18 

6.50 

0.050 

9.11 

7.79 

6.58 

5.85 

0.075 

9.92 

7.90 

5.63 

5.23 

0.100 

10.55 

7.96 

5.02 

4.72 

0.125 

9.21 

8.27 

4.84 

4.55 

0.150 

8.02 

8.38 

4.71 

4.39 

0.175 

7.22 

8.12 

4.59 

4.25 

0.200 

6.69 

6.92 

4.47 

4.14 
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FIG. 2. Electron-phonon coupling and critical temperature at 200 GPa of hydrogen chalcogenides with chemical substitutions. Upper 
panels: the integral function A'(cu). Lower panels: the electron-phonon coupling constant A (solid lines) and the critical temperature Tc 
(dashed lines) versus the substitution concentration x. The left, middle, and right panels respectively show the results for the case of VI-, V-, 
and Vll-substitutions. 


ical precompression effect^*^’^^’^^ on the atomic radius of VI 
elements. 

Further analyzing phonon linewidths reveals that they are 
larger for the H vibrational modes, particularly for the bond¬ 
stretching modes at high frequency, similar to the results by 
Errea et al}‘^ We note that the small magnitude of phonon 
linewidth in H 3 S 0 . 4 O 0.6 (Fig- Id) is ascribed to the reduc¬ 
tion of DOS. In order to expose more clearly the relative con¬ 
tributions of different phonon modes to the electron-phonon 
coupling constant A, we define an integral function A'(a;) = 
2 u}'~^a^F{ijj')dLo', where A'(oo) is A and a^F(w) is 
the Eliashberg function^^^®. As shown in Eig. 2a, the con¬ 
tribution from the intermediate frequency region (300 ^ 
1700 cm“^) has an increase in H 3 So. 7 Teo. 3 , which results 
in a corresponding increase in the electron-phonon coupling 
constant. Such a variation is a consequence of the increase 
of phonon linewidths in the intermediate frequency region, 
which is present in Eig. le but absent in the O- and Se- 
substitution systems (Eigs. Ic-ld). 

With the above results, we use equation (1) to estimate the 
Tc for the investigated Vl-substitution systems and the results 
are plotted in Eig. 2b. The effective Coulomb repulsion has 
been chosen to be p* = 0 . 12 , in the reasonable range of 
0.1 ~ 0.15^^, such that the estimated Tc for H 3 S and HsSe, 
194 K and 160 K, are most close to the reported values'®’^'. 
In all of the Vl-substitution cases, there is no substantial en¬ 
hancement of Tc compared with the 194 K in H 3 S, except a 
slightly higher Tc of 198 K in H 3 S 0 . 9 O 0 . 1 . Notably, despite 
a higher A in 11380 , 7 X 60 . 3 , the Tc is in fact reduced to 161 K. 
Because of the heavy mass of Te, the logarithmically aver¬ 
aged phonon frequency wm is only 90.8 meV in 11380 . 7 X 60 . 3 , 
much lower than the 125.2 meV in H 3 S. In fact, the effect of 


Vl-substitution at a fixed physical pressure can be viewed as a 
chemical pressure effect. Thus, our results that the Tc would 
decrease in the various Vl-substitution cases at 200 GPa are in 
agreement with the experimental observation'® and an earlier 
theoretical study^*, which have conhrmed that the Tc of H 38 
is peaked near 200 GPa. 

Besides the substitution of 8 by elements in the same main 
group, we also further consider substitutions by the adja¬ 
cent group elements, i.e., phosphorus and halogen groups. 
We mainly focus on the cases in which a small percentage 
(a; = 0.0 ^ 0.2) of S is replaced by P (Cl), and likewise Se 
replaced by As (Br). In sharp contrast to the Vl-substitution 
cases, the averaged phonon frequency is not expected to have 
a signihcant drop after substitutions, because the atomic mass 
of after chemical substitutions are close to that of the proto¬ 
types. Thus, as the calculations will show below, the changes 
of Tc follow closely the changes of A in the V- and VII- 
substitution cases. Note that the absence of imaginary fre¬ 
quency in the phonon spectra, as exhibited in Pig. 3 (also see 
supplementary materials), ensures the lattice dynamic stabili¬ 
ties for all the cases of interest. 

Due to the decrease of valence electrons after V- 
substitutions, the electronic energy bands shift up a little bit 
with respect to the Permi energy (see supplementary materi¬ 
als). Such changes can also be reflected effectively by the 
changes in DOS at the Permi surface, as summarized in Ta¬ 
ble I. It turns out that the DOS increase first and then decrease 
with increasing the V-substitution rates. Notably, the DOS can 
reach very large values, e.g., 10.55 in H3S0.9P0.1 and 8.38 
in H 3 S 60 . 85 As 0 . 15 , in units of Hartree~' per spin. Although 
the phonon spectra hardly change after V-substitutions, the 
phonon linewidths follow a close trend of the changes in DOS. 
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FIG. 3. Phonon spectra with phonon linewidth of V- and Vll-substitution systems, (a) H3S0.9P0.1. (b) H3So.gPo.2, and (c) H3So.gClo.2. 
The magnitude of phonon linewidths in H3S0.9P0.1 (H3So.gClo.2 ) is too large (small) and thus plotted with half (twice) of the real values. 


The magnitudes of phonon linewidths in the intermediate and 
high frequency regions rise sharply as the P-substitution rate 
reaches x = 0.1 and then fall steeply as the rate further 
increase to a; = 0.2, as shown in Figs. 3a-3b. Thus, the 
electron-phonon coupling constant A follows the same trend 
of the changes in DOS and in phonon linewidths, as shown in 
Fig. 2d. Analysis of A'(a;) plotted in Fig. 2c further confirms 
the enhancement of electron-phonon couplings in the interme¬ 
diate and high frequency regions in HsSo.gPo.i. Very similar 
behaviors can be found as the As-substitution rate increases 
from zero to x = 0.2 in H 3 Sei_xAsx, as seen in Figs. 2c-2d 
(also see supplementary materials). We find that the maximal 
coupling constant values are 2.12 in FlaSo.gPo.i and 1.97 in 
H3Seo.85Aso.15. 

Based on the above results, we use equation (1) with p* = 
0.12 to estimate Tc for the two V-substitution systems and plot 
them in Fig. 2d. The are found to be as high as 250 K in 
HsSo.gPo.i and 185 K in Fl3Seo.85Aso.15, which are greatly 
enhanced from the 194 K (160 K) in the prototypical FI3S 
(HoSe). Given the possible variation of the effective Coulomb 
repulsion p* in the range of 0.1 ~ 0.15, the T^. in HoSg.gPo.i 
may also vary in the range of 227 ~ 265 K. 

We now take into account the influence of varying the 
high pressure on the superconductivity*^’'® in the case of 
HoSi-xPx- Since we are particularly interested in optimiz¬ 
ing the Tc, we focus on two cases suggested by Fig. 2d: 
H3So.gPo.i and H3S0.g25P0.075- The main results are sum¬ 
marized in Table II and Fig. 4. In the case of H 3 S 0 .gP 0 .i 5 the 
Tc reaches the highest value at 200 GPa and remains the same 
value up to 225 GPa, whereas in the case of H3S0.925P0.0755 
the Tc increases monotonously as the pressure rises from 
150 GPa to 250 GPa, beyond which a structure instability is 
found. Results in the latter case follow the fact that the DOS 
at the Fermi energy, the phonon linewidths, and the electron- 
phonon coupling constants all increase gradually with increas¬ 
ing the pressure. For /r* = 0 . 1 , the Tc of H3S0.925P0.075 at 


250 GPa can reach 280 K, which is higher than the ice point. 
We note that the pressure of 250 GPa has already been fea¬ 
sible in the H 3 S experiment by Drozdov et a/*®. Given the 
range of fj,* = 0.1 ~ 0.15, the Tc in such an optimized case 
is at least as high as 241 K, as indicated by the error bar in 
Fig. 4c. (For n* = 0.12 that yields the Tc of 194 K in H 3 S, 
we find = 264 K.) 

Finally, in the Vll-substitution cases we have studied two 
cases: HoSi-xClx and HsSei-xBrx. In general, the elec¬ 
tronic energy bands shift down with respect to the Fermi 
energy, because of the increase of valence electrons; the 
DOS at the Fermi energy, the phonon linewidths, and 
hence the electron-phonon coupling constants all decrease 
monotonously as the degree of Vll-substitutions increases. As 
a consequence, the Tc drops monotonously with increasing 
the Vll-substitution degree. These results are summarized in 
Table I and Figs. 2e-2f (also see supplementary materials). 


III. DISCUSSION 

We have investigated the influence of partial atomic sub¬ 
stitution on the superconductivity of hydrogen chalcogenides 


TABLE II. Fermi-surface DOS (in units of Hartree”^/spin) and 
electron-phonon coupling constant A of HaSi-xPx under dif¬ 
ferent pressures. 


Pressure 

H 3 S 0 . 925 P 0.075 

H 3 S 0 . 

9 P 0.1 

DOS 

T ■ 



150 GPa 

9.16 

1.96 

9.13 

1.9 

175 GPa 

9.64 

2.0 

9.77 

1.98 

200 GPa 

9.92 

2.1 

10.55 

2.12 

225 GPa 

10.67 

2.24 

10.4 

2.19 

250 GPa 

11.1 

2.44 

9.87 

2.24 
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FIG. 4 . The effect of pressure in HaSi-xPx. (a)-(b) The phonon spectra and A'(tj) of H3S0.925P0.075 and H3S0.9P0.1 at 250 GPa. 
The magnitude of phonon linewidth is plotted with half of the real values. The solid (dash) lines plot A'(a;) at 250 ( 200 ) GPa. (c) Tc 
of H3S0.925P0.075 (red solid line) and H3S0.9P0.1 (black dash line) versus pressure. The error bars indicate the value range of with 
H * = 0.1 ~ 0 . 15 . 


using first-principles calculations with VGA. Our study has 
highlighted the key roles of strong covalent metallicity and 
low atomic mass in boosting the Tc of BCS superconductiv¬ 
ity. The former can produce large electron-phonon couplings, 
whereas the latter can yield high-frequency phonon modes, 
and in fact the highly-compressed H 3 S constitutes both ad¬ 
vantages. We now take H 3 S as the example to summarize our 
results. In the VTsubstitution cases, the reduction of DOS at 
Fermi energy dilutes the covalent metallicity even though the 
coupling constant can be enhanced in the Te-substitution case, 
the Tc is lowered because of the stronger atomic mass. This 
leads us to further study the cases of Cl- and P-substitutions, 
in which the atomic mass remains hardly changed. In the P- 
substitution case, the DOS, the phonon linewidths, the cou¬ 
pling constant, and hence the all increase as the substitu¬ 
tion rate increases from zero up to x = 0.1. In sharp contrast, 
the oppose trend occurs in the Cl-substitution case. 

In particular, we have shown that in the optimized case of 
H3S0.925P0.075 the Tc may reach a record high value of 280 K 
at 250 GPa, which is a feasible pressure in current experiments 
and would not induce any structure instability. In light of our 
results, it might be possible that the silicon-substitution would 
also enhance the Tc of high-pressure H3S superconductor. Be¬ 
cause of the even less valance electrons, the optimum Tc of 
HoSi-xSix might occur at even lower substitution concentra¬ 
tion, which might be easier to realize in experiment. For ex¬ 


ample, we hnd Tc = 274 K in the case of H 3 S 0 . 96 Si 0.04 (see 
supplementary materials). In the future, inclusion of anhar- 
monic effects'® may improve our predictions and thus better 
guide the experiments. Nevertheless, our hnding is exciting. 
It not only suggests that partial atomic substitution may lead 
to possible superconductivity above the ice point in the highly 
compressed H 3 S, but also gives a hope that in principle low 
atomic mass and strong covalent metallicity may be designed 
in novel materials to realize high-Tc BCS superconductivity 
under ambient pressures. 

IV. METHODS 

There are mainly two different procedures for disordered 
systems and partial atomic substitutions in hrst-principles 
calculations, i.e., an ordered supercell and the virtual crys¬ 
tal approximation (VCA)"'^. The former method is time- 
consuming and technically difficult to deal with the case for 
small concentrations. Thus, we chose to use VGA in this 
work. In calculations based on VGA, the primitive period¬ 
icity is retained but composed of virtual atomic potentials 
interpolating between the behaviors of actual components. 
Even though this approach neglects the local deformations 
around atoms and cannot explore the disordered structures 
very accurately, it often produces acceptable and useful re- 

































6 


suits that have been verified in many research fields of con¬ 
densed matter physics"^^’"^^^’^®’The atomic substitution 
in the present work is simulated by the self-consistent VCA. 
For example, in HsSi-^^Sej; the virtual pseudopotentials of 
the Si_a;Sea; is represented by a pseudopotential operator 
VvcA = a;Vse + (1 — a:)Vsj where Vse (Vs) is the pseudopo¬ 
tential of Se (S) atom. 

The present studies, including the electronic structures, the 
phonon spectra, and the electron-phonon couplings, were car¬ 
ried out using the ABINIT package®"^®^ with the local-density 
approximation (LDA). Hartwigsen-Goedecker-Hutter (HGH) 
pseudopotentials®^ were used in order to include spin-orbit 
couplings (SOC) for the heavy element, tellurium. The SOC 
for other elements were neglected since they are sufficiently 
light. By requiring convergence of results, the kinetic en¬ 
ergy cutoff of 30 Hartree and the Monkhorst-Pack fc-mesh 
of 40x40x40 were used in all calculations about the elec¬ 
tronic ground-state properties. The phonon spectra and the 
electron-phonon couplings were calculated on a 8x8x8 q- 
grid. Since calculating electron-phonon couplings referred to 
integrals over the Brillouin zone, we also carefully checked 
convergence of the results on the aforementioned fc-mesh and 
g-grid, by comparing them with results in denser samples (a 


40x40x40 fc-mesh and a lOx lOx 10 g-grid). 
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